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ABSTRACT 


A  review  of  the  current  theory  and  methods  for  the  computa- 
tion of  maximum  flow  in  networks  is  presented  along  with  a  simplified 
method  for  determination  of  a  feasible  flow  in  networks  with  upper 
and  lower  bounded  arcs.  A  computational  procedure  is  presented  which 
is  used  to  calculate  the  maximum  flow  for  a  general  network.  The 
network  is  reduced  to  an  equivalent  basic  network.  An  associated 
network  is  used  to  compute  a  feasible,  then  the  maximum  flow  for 
the  basic  network.  A  computer  program  is  included  for  use  in 
computation  of  maximal  flows  in  large  networks. 
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CHAPTER  I 

INTRODUCTION 

The  growing  uses  for  mathematical  programming  models  1n  trans- 
portation and  communication  networks  have  prompted  study  in  the  devel- 
opment of  computational  methods  to  determine  maximal  network  flows. 
The  subject  of  this  thesis  is  the  mathematical  problem  of  determining 
maximal  steady  state  flows  in  networks  which  are  subject  to  capacity 
limitations  on  the  arcs  and  nodes  of  the  network. 

An  efficient  and  widely  used  method  to  compute  maximal  flow 
was  developed  by  Ford  and  Fulkerson  [7].  This  method  along  with  a 
review  of  the  current  theory  will  be  presented  in  Chapter  II  in 
order  to  give  the  reader  the  theoretical  foundation  upon  which  the 
remainder  of  the  thesis  is  based.  Chapter  III  introduces  a  new 
method  which  utilizes  the  Ford  and  Fulkerson  algorithm  to  compute  the 
maximal  flow  in  networks  with  non-zero  lower  bound  limitations  on  the 
flow  in  the  arcs.  Since  many  networks  of  interest  are  more  general 
in  nature,  the  generalized  network  is  introduced  in  Chapter  IV  along 
with  a  computational  procedure  for  the  determination  of  maximal  flow. 
To  illustrate  the  techniques  proposed,  an  example  is  included.  The 
procedure  can  be  used  to  solve  small  problems  by  hand,  but  the  amount 
of  work  increases  rapidly  with  the  size  of  the  network.  Therefore,  a 
computer  program  which  may  be  used  to  solve  larger  problems  is 
included  in  the  appendices.  A  brief  description  of  the  program  is 
found  in  Chapter  IV. 


Definitions  and  Symbols 

In  order  to  establish  a  common  understanding  as  to  the  exact 
meaning  of  various  terms  and  symbols  used  to  present  the  material  of 
this  thesis,  an  initial  set  of  definitions  is  presented.  Other  terms 
will  be  defined  throughout  the  material  as  needed  to  facilitate  the 
presentation. 

Define  a  network,  6,  as  a  graph  which  can  be  represented  on 
a  plane  in  such  a  way  that  the  set,  N,  of  nodes  are  distinct  points 
and  the  set,  A,  of  arcs  are  simple  curves  which  connect  two  distinct 
nodes.  Furthermore,  no  two  arcs  can  meet  except  at  the  nodes  which 
are  their  extremities.  The  graph  contains  no  loops. 

Define  uT(N)  to  be  the  set  of  arcs  or  flows  incident  to  a  set 
of  nodes,  N,  and  eo  (N)  to  be  the  set  of  arcs  or  flows  incident  from  N. 

A  directed  arc,  a.  .,  is  defined  as  an  arc  incident  from  a  node 
i  and  incident  to  a  node  j. 

An  undirected  arc  connects  node  i  and  node  j  without  having  an 
orientation. 

Flow  can  be  thought  of  as  a  value  of  the  steady  state  rate  of 

movement  of  a  homogeneous  commodity  along  a  path  or  channel.  Positive 

flow  in  a  directed  arc,  a.  , ,  will  move  with  the  orientation  of  the 

arc  and  will  be  denoted  by  the  symbol,  y,  ...  Positive  flow  in  an 

undirected  arc,  denoted  y.  .,  may  move  in  either  direction  but  not 

in  opposite  directions  at  the  same  time. 

To  each  arc,  a.  .  c  A,  there  will  be  associated  two  real 
•  »J 

numbers,  b.  .,  and  c,  ,,  which  represent,  respectfully,  the  minimum 
•  »j      i  »j 

and  maximum  allowable  flow  1n  that  arc.  The  bounds  associated  with 
undirected  arcs,  a..  .,  are  denoted  b\  .,  and  c.,  .. 
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An  arc  is  saturated  if  y^  .  ■  c^  ... 

Define  Y_  as  the  set  of  flows  in  a  network  G. 

Initially,  the  capacities  of  the  nodes  will  be  assumed  to  be 
infinite.  The  nodes  will  be  classified  according  to  the  value  of  the 
difference  of  the  sum  of  the  flows  Into  the  node  minus  the  value  of 
the  sum  of  the  flows  out  of  the  node.  A  node  1s  designated  as  a 
source  node  1f  this  difference  1s  negative,  and  as  a  sink  node  1f  this 
difference  1s  positive.  For  all  Intermediary  nodes,  the  sum  of  the 
flows  into  a  node  must  equal  the  sum  of  the  flows  out.  The  set  of 
source,  Intermediary  and  sink  nodes  will  be  denoted  as  S,  I,  and  T 
respectfully. 

Define  F  as  the  value  of  the  flow  1n  a  network. 

A  flow  is  feasible  1f  and  only  1f: 

(1)  bi j ± yi ,j ± ci ,j     Vai.j£A 

(2)     F  -  I  4      y, 

ys.ji 


,c/(S)  S'J 


E  .      y,  t  -  F  =0 

y1ftcu)(T)  1'< 

(3)  S  .      y,  1  -  E  yi   1  o  o        1,j  c  I 

Define  a  chain  as  a  sequence  of  arcs  such  that  each  arc,  a.  . 
ij  c   I,  1s  connected  to  an  arc  a.  ,  at  node  1  and  to  arc  a.  ,  at  node  j 
A  cycle  1s  a  sequence  of  arcs  where: 

(1)  The  sequence  1s  a  chain. 

(2)  The  sequence  does  not  use  the  same  arc  twice. 

(3)  The  Initial  and  terminal  nodes  of  the  chain  coincide. 


A  cut  1s  defined  as  any  set  of  directed  arcs  containing  at 
least  one  arc  from  every  chain  of  positive  capacity  joining  the 
source  to  the  sink.  The  value  of  the  cut  1s  equal  to  the  sum  of  the 
capacities  of  the  arcs  of  the  cuts. 

Basic  Network 

Much  of  the  Initial  development  of  the  material  1n  this  thesis 
deals  with  computation  of  flows  1n  simple  or  basic  networks.  There- 
fore, for  the  purpose  of  this  paper,  a  basic  network,  G,  1s  defined 
as  a  network  which  has  the  following  properties. 

1.  The  network  has  a  return  arc,  an  •,,  which  1s  added  to  the 
network  for  computational  purposes  only. 

2.  Node  1  has  only  one  arc,  an  i »  Incident  to  1t.  Node  1  1s 

n,  i 

designated  as  the  source.  In  the  basic  network  there  1s 
only  one  source  node. 

3.  Node  n  has  only  one  arc,  a_  ,,  Incident  from  1t.  Node  n 

n,  i 

1s  designated  as  the  sink.  In  the  basic  network  there  1s 
only  one  sink  node. 

4.  For  each  arc,  a^  j,  the  associated  c,  ,,  1s  positive  and 

bi.j  ■  °- 

5.  Two  distinct  nodes  (1,J)  may  be  connected  by  only  one 
directed  arc.  The  network  contains  no  undirected  arcs. 
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CHAPTER  II 

REVIEW  OF  THE  LITERATURE 
In  order  to  acquaint  the  reader  with  the  subject  of  maximal 
flow  in  networks,  a  review  of  the  current  theory  and  primary  combina- 
tional techniques  is  presented  in  this  chapter. 

Basic  Network  Theory 

A  great  deal  of  theoretical  work  in  networks  has  already  been 
done.  A  few  of  the  relevant  theorems  are  presented  here  in  order  to 
provide  justification  and  a  rationale  for  the  combinatorial  techniques 
and  algorithms  that  are  presented  in  later  sections. 

Theorem  1.  Opposite  directed  flows  on  the  same  arc  cancel  [4] 

This  theorem  means  that  given  an  undirected  arc,  a,  .,  with  flows 

•  »J 

y.  .  >  0  and  y.  .  >^  0,  these  two  flows  can  be  replaced  by  the  flows 

y'i,j and  y,j.i 

Theorem  2.  A  set  Y  with  F>0,  satisfying  the  capacity  con- 
straints and  node  conservation  equations  can  be  decomposed  into  a  set 
of  positive  chain  flows  from  the  source  to  the  sink  and  a  set  of 
circular  flows  such  that  the  direction  of  positive  flows  in  any 
common  arc  is  the  same  for  all  chains  [4]. 

Theorem  3.  There  exists  a  positive  flow  from  the  source  to 
the  sink  if  there  exists  at  least  one  chain  of  arcs  with  positive 
capacity  from  the  source  to  the  sink  [4]. 

Theorem  4.  A  flow  F°  is  maximal  if  and  only  if  the  maximal 
flow  is  zero  in  a  second  network  formed  by  replacing  c.  .  by  c  .  - 
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Theorem  5.  Given  any  partition  of  nodes  into  two  sets,  where 

the  first  set  includes  the  source  node  and  the  second  set  the  sink, 

then  a  feasible  solution  F  is  maximal,  if  every  arc,  a.  .,  that  joins 

1  »j 

a  node  in  the  first  set  to  a  node  of  the  second  is  saturated  [4]. 
Theorem  6.  The  maximum  flow  in  a  network  is  equal  to  the 
minimum  cut. 

Tree  Method 

Dantzig  [4]  has  developed  a  systematic  procedure  for  com- 
puting the  maximal  flow  in  a  network  by  using  chains  connecting  the 
source  and  the  sink  of  the  network.  A  knowledge  of  the  tree  method 
will  give  a  better  understanding  of  the  Ford  and  Fulkerson  algorithm 
which  is  presented  in  the  next  section  because  the  two  methods  are 
identical  in  principle. 

A  tree  [1]  is  defined  as  a  connected  graph  with  at  least  two 
nodes  and  no  cycles. 

Consider  a  basic  network  such  that: 

bifj  =  0,cu>0      Va1>j£A 

The  arcs  may  be  directed  or  undirected.  Initially  all  arcs 
are  unmarked  and  the  feasible  flow  is  zero.  In  order  to  determine 
the  maximum  flow: 

1.  Create  a  tree  such  that  the  arcs  are  unsaturated. 

2.  Select  two  sub  trees  -  one  T, ,  branching  out  from  the  source, 
node  1,  and  the  other,  T  ,  branching  out  from  the  sink,  n,  such  that 
each  node  1s  reached  by  just  one  arc  of  the  tree. 

3.  Choose  one  unsaturated  arc  which  connects  the  two  trees.  Thus, 
there  will  be  just  one  chain  from  node  1  to  node  n.  If  none  can  be 
found,  go  to  step  7. 
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4.  Determine  the  amount  that  the  flow  can  be  increased  along  the 
chain  from  node  1  to  node  n,  by  the  following  procedure.  Determine 
the  amount  that  the  flow  can  be  increased,  AX,  where  AX  =  min  (c  ,  - 
y-  •»  y^  i)'  Flow  along  the  chain  directed  toward  the  sink  is 
designated  y.  .  and  flow  directed  toward  the  source  along  the  chain 
is  designated  y.  , .  There  are  two  cases: 

(a)  AX  >  0,  increase  the  flow  by  AX  along  the  chain,  go  to 
next  step. 

(b)  AX  =  0.  This  means  one  arc  in  the  chain  is  saturated 
by  a  flow  directed  toward  the  sink,  go  to  next  step. 

5.  Mark  one  of  the  saturated  arcs. 

6.  Eliminate  the  saturated  arc  from  the  chain.  Thus,  there  are  two 
sub-trees.  Return  to  step  3. 

7.  Check  all  marked  arcs  joining  T,  to  T  .  If  all  these  arcs  have 

flow  from  T,  to  T  then  the  optimum  has  been  reached.  If  any  arc 

has  flow  from  T  to  T, ,  use  this  arc  to  connect  T,  to  T  and  proceed 
n    I  I    n    r 

with  step  4. 

Ford  and  Fulkerson  Algorithm 

Given  the  basic  network,  the  most  direct  method  developed  to 
date  of  computing  the  maximum  flow  is  the  algorithm  developed  by  Ford 
and  Fulkerson.  This  method  is  the  one  most  cited  in  the  literature 
concerning  maximal  network  flow.  The  procedure  is  outlined  below. 

1.  Begin  with  a  feasible  flow  in  0  in  all  arcs. 

2.  Label  each  node  as  follows: 

a.  Label  node  1,  with  the  number  1. 
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b.  If  node  i  is  labeled  and  node  j  is  not  labeled,  label 

node  j  with  the  number  (i)  if: 

(1)  a.  ,  exists  and  y.  .  <  c.  . 

I  »j  i  »j    i  »j 

(2)  a .  .  exists  and  y .  .  >  0 

c.  If  the  sink,  node  n,  becomes  labeled  by  this  procedure, 
then  the  flow  from  node  1  to  node  n  can  be  increased. 

3.  If  node  n  is  labeled,  construct  a  simple  chain  from  node  1 
to  node  n  by  back  tracking  from  n  to  1  using  the  labels 

on  the  nodes. 

a.  The  chain  created  will  have  the  property  that  all  arcs, 

a.  .  on  the  chain  directed  toward  the  sink  have  a  flow 
■  »J 

0  <  y.  .  <  c.  .  and  all  arcs  a,  .  directed  toward  the 
source  will  have  a  flow  y,  ;  >  0. 

b.  The  flow  can  be  increased  by  an  amount  AX  where: 
AX«m1^c.ja.  -yui,  yljk)   V  »uy   a,  >k  in  chain. 

4.  Increase  the  flow  along  the  chain  by  an  amount  AX  and 
return  to  step  2. 

5.  This  procedure  is  repeated  until  the  sink  cannot  be 
labeled.  Flow  is  maximized  when  this  is  the  case. 

Now  consider  the  more  general  problem  in  which  the  lower  bounds 

on  the  directed  arcs  are  no  longer  assumed  to  be  zero.  Thus  the  flow 

in  each  arc  is  bounded  both  above  and  below. 

Such  that:  -  °°  <_  b.  .  <_  y .  .  <_c.    .   <_  +  °°    a.  .  £  A 
Theorem  1  ■;     In  the  case  of  networks  with  lower  bounds  on  the 

arcs,  C  is  defined  as  the  subset  of  N  such  that  if  node  1  k   C  the  node 

n  ^  C  or  if  node  1  c_  C  then  node  n  £  C. 
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A  flow  y.  .  exists  and  is  feasible  if  and  only  if: 
i  >J 

!-  b1,j^1.j±c1.J   Vai,j^A 

2.      Z      c.  ,  >     I  b.  , 

a.  .  c  W-(C)  1jJ   a.  .  c  </(C)   1  'J 
i  »j  •  »j 

This  more  general  type  of  problem  can  be  satisfactorily  worked 
by  using  the  Ford  and  Fulkerson  algorithm  if  a  feasible  flow  can  be 
found. 

The  Ford  and  Fulkerson  algorithm  can  be  used  to  solve  for  a 
maximum  flow  in  an  associated  network  G'  which  is  derived  from  the 
original  network,  G.  This  method  is  outlined  below. 

1.  Create  an  associated  network  G'  such  that  G'  contains 

the  nodes  and  arcs  of  the  original  network  G  with  arc 

capacities  c'.  .  =  c.  .  -  b.  .  .  The  lower  bound, 
i  >J    >  »J    i  >J 

b'.  j  equals  zero. 

2.  Add  a  new  source,  labeled  0,  and  a  new  sink,  labeled  m, 

(m  =  n  +  1),  to  G1. 

a.  If  the  arc,  a.  . ,  exists,  then  construct  two  new 
i  >J 

arcs  according  to  the  value  of  the  lower  bound,  b.  .. 

If  b4  .  >  0,  then  construct  the  arcs,  a'.   of  capacity 
i  ,j  l  ,m 

c'j  m  =  b.  .  ana  a'   .  with  capacity  c'   ,  =  b.  .  . 
i  ,m    i  ,j      o,j       K   J   o,j    i , j 

If  b.  .  <  0,  then  construct  the  arcs,  a'   i  >/ith 
i  »J  0,1 

capacity  c*  .  =  -  b.  .  and  a'.   with  capacity 
o,i      i ,  j      j  ,m 

c'  .   =  -  b.  .  . 
J,m     l  ,j 

bo  The  Ford  and  Fulkerson  algorithm  can  then  be  used  to 
determine  the  maximum  flow  in  G'  from  node  0  to  node  m. 

3.  Two  cases  can  result  from  step  2b. 

a.  Case  one:  If  the  flow  in  each  arc,  a'   .  and  a' , 

o,j      j,m 
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j  c_  N',  equals  the  capacity  of  each  arc,  then  the  flow 

is  feasible  for  the  network  G  and  the  problem  has  a 

solution. 

b.  Case  two:  If  for  any  arc,  a'   .  or  a' .   j  c  A' , 

o, j      j  ,rn 

the  flow  in  that  arc  does  not  equal  the  capacity  of 
the  arc,  the  problem  has  no  solution. 

4.  If  a  feasible  flow  for  G  can  be  found  by  the  above  pro- 
cedure, transfer  the  flow  from  G'  to  G  by  using  the 
following  transformation: 

y.  .  =  y' .  .  +  b.  .   V  a.  .  c  A 

5.  Ford  and  Fulkerson's  algorithm  is  then  used  to  determine 
the  maximum  flow  using  the  feasible  flow  as  a  starting 
basis. 
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CHAPTER  III 

A  SIMPLIFIED  METHOD  FOR  COMPUTATION  OF  MAXIMAL  FLOW 
To  compute  a  maximum  flow  in  networks  with  non-zero  lower 
bounded  arcs,  an  initial  feasible  flow  must  be  determined.  The  method 
presented  in  the  last  chapter  produces  a  complex  associated  network, 
G',  even  with  a  relatively  simple  network,  G.  The  principle  advantage 
of  the  method  presented  here  [10]  is  that  an  initial  feasible  flow 
can  be  determined  by  use  of  a  more  simple  network,  G',  than  in  the 
previous  method. 

Proposed  Algorithm 

The  development  of  the  algorithm  is  based  upon  a  consequence 
of  the  flow  balance  equations.  In  order  to  provide  a  rationale  for 
the  procedures  presented,  the  development  of  the  associated  network 
will  be  presented  now. 

Given  a  basic  network  G  where  the  flow  in  each  arc  is  bounded 
as  follows: 

-»±bi,jiyi>jicilj+"    YjCA 
The  problem  of  determining  the  maximum  flow  in  G  can  be  stated  as 
follows: 

Maximize:  y  , 

Subject  to:      I  y.  .  -     Z     y.  .  =  0  V  i  c  N 

y,  ,  cuf(i)  J'1   y.  .  cc/(i)  !-J 
j » '  ■  »j 

bi,j±yi,j±cij     Vau^A 

Now  transform  the  entire  network  G  into  an  associated  network 
G'  using  the  change  of  variables: 

yi,o=y'i,j  +  bi,j   Vai,j^A 
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the  problem  then  becomes: 

Maximize:  y'n  1 

Subject  to: 

E     y' .  .  -      I  y'-  •  +     E     b.  j 

ID  J>1   y'-,co)(i)  1,J   b.  .  c  uf(i)  J'1 


^.J*"*1*      ^.j^^       Dj,i 

E     b.  .  =  0    V  i  c  N 

b,  ,  cu3+(i)  lsJ 
1  »j 

and:  0  <  y' .  .  <  c1 .  .  =  C  .  -  b.  • 

Let     E     b.  .  -     E     b.  .  =  d.     i  c  N 
b.  .  c  a)-(i)  Jsl   b.  .  co)+(i)  lsJ    n 

Note  two  things.  First  the  sum,   E  d. ,  is  equal  to  zero. 

ic  N  n 
Second,  if  a  flow  of  -d.  was  passed  through  each  node,  the  problem 

would  be  exactly  the  problem  of  flow  maximization  in  a  basic  network. 

Hence,  a  technique  similar  to  Berge's  [2]  can  be  used  in  which  new 

nodes  0  and  m  =  n+1  are  added  to  the  network  G*.  Therefore,  only  a 

maximum  of  one  arc  per  node  is  added  to  G'  in  this  method. 


Procedure 

Using  the  results  of  the  previous  section  the  solution  to  the 
maximal  flow  in  the  above  network  can  be  determined  using  the  following 
procedure. 

1.  Add  a  new  source  and  sink  node,  denoted  0  and  m  (m  =  n  +  1 ) , 

S 
respectfully  to  G' . 

2.  The  arcs  a'   .  and  a1 .  m  are  added  as  follows.  If  d.  >  0, 

o,i  j,m  l         ' 

tl^an  arc,  a'      .,  is  constructed  with  capacity  d..     If  d.  <  0,  an 

0,1 

arc,  a'.  ,  is  constructed  with  capacity  -d. .  If  d.  -  0,  no  arc  1s 
added  to  G'  for  node  i. 
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3,  The  Ford  and  Fulkerson  algorithm  is  then  used  to  maximize 
flow  in  G'  from  node  0  to  node  m.  Again  as  in  the  previous  method 
this  maximum  flow  is  equivalent  to  a  feasible  flow  in  the  network  G 
if  and  only  if: 

y'o>.  .  d.  V13d>0 

y'j,m  =  di  Vi3d<0 

4,  If  the  problem  is  solvable  continue  by  eliminating  all 
arcs  incident  from  node  0  and  incident  to  node  m.  Node  1  and  node  n 
are  now  the  source  and  sink  respectfully  for  G'. 

5,  Using  the  flow  generated  in  step  3  as  a  starting  basis, 
maximize  flow  from  node  n  to  node  1  by  again  using  the  Ford  and 
Fulkerson  Algorithm. 

6,  Upon  completion  of  step  5,  transfer  flow  from  G'  to  G  by 
use  of  the  transformation  equation: 

v.  .  ■  y1 .  . ;  +  b.  .    V  a.  .  c  A 

The  maximum  flow  in  G  is  equal  to  y  , . 

n  >  i 
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CHAPTER  IV 

PROCEDURE  FOR  THE  SOLUTION  TO  THE  GENERALIZED  MAXIMUM  FLOW  PROBLEM 

A  proposed  methology  for  the  computation  of  a  solution  to  the 
generalized  maximal  flow  type  problem  as  defined  below  is  presented. 
The  generalized  network  is  first  reduced  to  a  basic  network  and  then 
the  technique  presented  in  Chapter  III  is  used  to  solve  for  the  equiv- 
alent maximum  flow  in  the  basic  network.  An  example  is  presented  in 
order  to  demonstrate  the  techniques  discussed  in  this  chapter.  A 
description  of  the  computer  program  contained  in  Appendix  A  for  use 
in  solving  maximal  flow  problems  for  large  networks  is  included. 

The  Generalized  Network 

Generally  maximal  flow  solutions  are  required  for  networks 
which  do  not  fit  the  restrictive  definition  of  a  basic  network.  For 
the  purpose  of  this  thesis  the  generalized  network  is  defined  as  a 
modified  basic  network  with  a  set  of  arcs  A,  and  a  set  of  nodes  N, 
such  that: 

1.  Each  set  uf(1)  and  u>  (i)  can  contain  any  number  of 
directed  and  undirected  arcs. 

2.  Let  D  be  a  subset  of  A  such  that  if  and  only  if  the  arc 
a.  .  is  directed  it  is  a  member  of  D. 

3.  Let  U  be  a  subset  of  A  such  that  if  and  only  if  the  arc 
a.  .  is  undirected  it  is  a  member  of  U. 

4.  K  is  a  subset  of  N  such  that  if  the  node  i  is  bounded 

it  is  a  member  of  K.  To  the  bounded  node,  i,  we  associate 
two  real  numbers,  b.,  the  lower  bound,  and,  c,  the  upper 


20 


bound  such  that  the  total  flow  through  the  node  i, 
denoted  y. ,  is  bounded: 

bi  i  yn- 1  ci       1  £  K 

5.  All  nodes,  j,  not  bounded  belong  to  the  subset  1,  where 

L  c  N. 

6.  The  network  may  or  may  not  contain  multiple  sources  or 
sinks,  but  shall  contain  at  least  one  source  and  one  sink. 

Procedure  for  Reduction  of  the  Generalized  Flow  Problem 

In  order  to  solve  for  the  maximal  flow  in  the  generalized  net- 
work by  the  techniques  presented  in  Chapter  III,  the  network  must  be 
transformed  into  an  equivalent  basic  network. 

The  proposed  methods  for  reduction  are  based  upon  results  of 
the  general  theorems  presented  in  Chapter  II.  References  cited  give 
further  proof  of  some  of  the  methods  used. 

First,  consider  the  case  in  which  the  generalized  network,  G, 
contains  both  directed  and  undirected  arcs.  This  network  can  be 
simplified  by  replacing  each  undirected  arc  by  two  directed  arcs  as 
follows: 

Given  an  undirected  arc  a-  .  £  U 

Replace  a.  A   by  a.  .  and  a.  .  c  D 

Where 

0  <  y,  .  <  c.  .  =  c.  . 

0  <  y .  .  <  c.  ,=c.  , 

Upon  reaching  a  solution  for  the  maximal  flow,  the  orientation 

by  the  flow  in  the  undirected  arcs  can  immediately  be  determined  from 

the  values  of  y.  4  and  y.  .  since  one  flow  will  be  zero. 
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Now  consider  a  network  in  which  two  distinct  nodes,  i,  j,  are 
connected  by  more  than  one  arc.  Successive  arcs,  flows,  and  bounds 
are  denoted  by  prime  symbols.  For  example  if  there  are  m  arcs  from 
node  i  to  node  j,  the  kth  arc  would  be  denoted  a.  /  ' .  It  will  be 
assumed  that  each  arc  has  both  upper  and  lower  bounds. 

b.^<y.(k)<c(k)       k  =  l,2,...m 

The  set  of  arcs  a.  .'  '  k  =  (l,...m)  may  be  replaced  by  a  single  arc 

m    (k)        k      m  k  (k) 
a.  ,,  where;  c.  .  =  E  cV        b,  .  =  -  v  b.  .' 

1  ,J  1  ,J     .   ,    ■  sj  I  »J     ^=1    I  »J 

Next  assume  a  pair  of  nodes  (i,j)  are  connected  by  two  oppositely 

directed  arcs  a.  .  and  a.  .  c  D. 
i  ,J     J.i  - 

The  pair  of  arcs  may  be  replaced  by  a  single  directed  arc 

in  either  direction  as  follows:  If  c.  .  >_  b.  . 

Replace  a.  .  and  a.  .  by  a.  .  such  that:* 
l » J     J  > '     i  >J 

c' .  .  =  c.  .  -  b.  . 
i.J    i.J    "i.J 

b' .  .  =  b.  .  -  c.  . 
i.J    i.J    t.J 

Otherwise  orient  the  arc  from  j  to  i. 

If  a  network  has  multiple  sources  and,  or  sinks,  merely  create 
a  single  source  and,  or  single  sink.  Then  construct  one  directed  arc 
of  infinite  capacity  from  the  new  source  to  each  old  source  and  one 
directed  arc  of  infinite  capacity  from  each  sink  to  the  new  sink  [8]. 

If  the  flow  through  the  nodes  is  restricted  to  lie  between 
two  values,  the  character  of  the  maximum  flow  problem  itself  remains 
unaltered.  Assume  we  are  given  a  network  in  which  the  flow  is  subject 
to  both  arc  and  node  capacities,  where  all  arcs  are  directed  or  have 
been  converted  to  directed  arcs  by  the  previous  techniques.  This 
network  can  be  reduced  to  the  basic  network  maximum  flow  problem  by 
the  following  procedure.  Expand  the  network  such  that  for  each  bounded 
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node  (1)  £  K,  two  n.  v  . abounded  nodes  (i1,  fM)  r?  created,  so  th?t 
the  single  directed  arc,  a.,  .,,  ,  joins  the  two  new  nodes  with  the 
upper  and  lower  boun  is  associated  with  the  bounded  node  in  the  reduced 
network.  To  each  arc,  a.  .,  create  an  arc  a'.  .,  in  G'  and  to  each 
arc  a.  .  in  G  create  an  arc  a.,,  .  in  G'.  The  bounds  on  the  other 
arcs  remain  the  same. 

Example  Using  the  Simplified  Technique 

Given  the  generalized  network  chowi  in  Mgure  1,  we  are 
required  to  compute  the  maximal  flow  from  node  2  and  3  to  node  10.  Or; 

Maximize  F-     l  J     y1>n  •  y5>10  +  y7(10  +  y6>1„ 

'1  ,n  — 

The  flow  in  each  arc  is  bounded  as  follows: 

3  1*2,4  I8  0<y256<4         3<y3>4<7 

5  <  y3  5  <  10  3  <  y4  5  <  8  0  <  y4  5  <  5 

(1)  (2) 

0iy4,6l3  ]  1*4,6  I2  °  1*4,6  I2 

3<y5>4<8  0<y5j6<8  3  <  y&_  ^  <  7 

2i*5,ioi5        ]  1*6,7  I3       31^6,lol9 

(1)  (2) 

1  1*7,5  I4  ]  1*7,6  I3         2  ly7,6  I5 

4<y7JO<  10 

The  flow  in  all  nodes  except  5  and  6  is  unrestricted.  The  flow  in 
node  5,6  is  restricted  as  follows, 

4 1  ys  i 15   4 1  ye  i  i° 

This  flow  problem  can  be  simplified  by  using  the  techniques 
presented  in  the  last  section.  The  procedure  is  outlined  by  ordered 
steps.  Although  many  steps  could  be  combined  for  this  simple 
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FIGURE  1 

GENERALIZED  NETWORK  G  WITH  UPPER 
AND  LOWER  BOUNDS  ON  ARCS  AND  NODES 


FIGURE  2 

GENERALIZED  NETWORK  WITH  SINGLE 
ARCS  REPLACING  MULTIPLE  ARCS 
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example,  it  is  felt  that  for  hand  and  computer  computation  this  is  the 
most  efficient  procedure. 

1.  Add  the  single  source  node,  denoted  as  1,  to  the  network. 
In  addition,  add  the  arcs  a,  2  and  ai  3  to  the  network 
such  that: 

0  1  y-\  2  -  +  °° 
0  <y1}3  <  +  «> 

2.  Replace  each  undirected  arc  a.  .  with  upper  bound  c.  .  >_  0 
by  two  directed  arcs  a.  .  and  a.  . . 

Therefore: 

(2)  (2)      (2)         (?) 
a4j5  becomes  a^£  and  a£^  0  <  y^£  <  5,  0  <  y£^  <  5 

(3)  (3) 

a4,6  becomes  a4 ,6  and  a6,4  °iy4,6^2'  °  1^6,4  ^2 

1         55,6  becomes  a5j6  and  a6>5  0  <  y5j6  <  8,  0  <  y6>5  <  8 

3.  Combine  all  arcs  joining  each  pair  of  arcs  (1,j)  according 
to  the  following  procedure. 

a.  Combine  all  arcs  which  are  similarly  directed,  into  a 
single  arc.  Thus  the  arcs  ai  I,   ai  i,  ai  I   become  a.  g 

c    =  c(1)  +  c(2)  +  c(3)  =3+2+2=7 
c4,6   C4,6   C4,6   C4,6   J  d       6       ' 

h   =  h'  '  +  vS  '   +  h'  '  =  i  +  n  +  n  =  1 
b4,6   D4,6   D4,6  +  D4,6   '   U   U   ' 

Similarly  form  the  arcs  a4  c>  an  4  and  a7  c      3 

c4,5  =  13  b4,5  =  3 
c5,4  =  13  b5,4  =  3 
c7,6=  8        b7,6=8 


9 
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FX6URE  3 
EQUIVALENT  BASIC  NETWORK 


FIGURE  4 

TRANSFORMED  NETWORK  G' 
26 


b.  Transform  each  pair  of  arcs  which  are  directed  in 

opposite  directions  between  each  pair  of  nodes 
into  a  single  arc. 

Thus:  a,  5  and  a5  *  become  a*  5  3 

c4,5  =  c4,5  "  b5,4  "W-3-10 

b4,5  "  b4,5  "  c5,4  =  3  -13  «  -  10 

Similarly  the  pair  of  arcs  a.  r   and  a,  ,;  a5  6  and 

a6  5*'  a5  7  anc*  a7  5;  a7  6  anc'  a6  7  are  combined.  Tne 
simplified  network  is  shown  in  Figure  2. 

3.  Complete  the  transformation  to  the  basic  network  by  the 
following  procedure. 

a.  Expand  the  network  so  that  all  nodes  have  infinite 
capacity  for  flow  as  outlined  earlier.  Thus  for 
nodes  5,  and  6  the  new  arcs  a5  g  and  ag  g  are 
created  and  the  equivalent  network  shown  in  Figure  3 
is  established. 

b.  Add  the  return  arc  a-|Q  ■,  for  computational  purposes. 

4.  Transform  the  network  G  to  G'  by  using  the  change  of 
variables  equation   y.  .  =  y'.  .  +  b.  .    a.  .  £  A 

*UT,1  =^10,1  ^1,2  =y'l,2  ^1,3  =  y'l,3 

*2,4  =^2,4  +8  *2,6  ■^2.6  y3,4  =  ^3,4  +  3 

*3S5  =  ^3,5  +  5  *4,5  =  ^4,5  "  10  *4,6  =  ^4,6  "  ] 

^5,8  =^5,8  +4  ^6,9  =^6,9  +  4  *7,6  =  ^7,6 

^7,10  =  y'7,10  +  4  *8,6  =^8,6  "  8  ^8,7  =  ^8,7  "  ] 

^8,10  =  ^'8,10  +  2  ^9,10  =  ^9,10+  3 
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The  equivalent  problem  is  obtained: 

Maximize  y',Q  ,  with  the  node  flow  satisfying  the  balance 

equations  as  follows: 


Node 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 


10,1 

1,2 

1,3 

3,4 

3,5 

2,6 

8,7 

5,8 

6,9 

8,10 


-  y 

-  y 

-  y 
+  y 
+  y 
+  y 

-  y 

-  y 

-  y 


1,3 

2,4 

3,4 

2,4 

4,5 

4,6 

7,6 

8,6 

9,10 

7,10 


-*'l,2 

"y,2,6 

-y'3,5 

-^4,5 

-^5,8 

+y,7,6 

-^7,10 

"  y,8,7  '  y,8,10 
+  1 


-  8 

-  8 

"y,4,6 

-  9 

+y,8,6 

-  5 


+  22 


-  y 


6,9 


-  13 


9,10 


-  y 


10,1 


+ 11 


+  9 


=  0 

=  0 
=  0 
=  0 
=  0 
=  0 
=  0 
=  0 
=  0 
=  0 


5. 


Therefore  the  values  of  the  remainder  d.  are  immediately 

seen  to  be  as  follows: 

d]  =  0  d2  =  -  8  d3  =  -  8  d4        =  +  22 

d5  =  -  9  d6  =  -  13  d?  =  -  5  dg  =  +11  dg  =  +1  d]()  =  +  9 

Using  the  values  of  d. ,  i  =  i,,,,10,  as  the  capacities,  the 

new  source,  sink  and  arcs  are  added  to  obtain  the  equivalent 

network  G'  shown  in  Figure  4. 

Then  using  the  Ford  Fulkerson  algorithm  to  obtain  a  feasible 

flow  in  G  the  following  values  are  obtained  for  the  flow  in 

G\ 
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y'l,2     =  12       y'l,3=8    ^2.4     '     °     V'2,6     =  4 

y'3>4   ■   o    y,§B-o  y4§s   -14   y4>6   =8 

y'5.8  =  5   y,6,9  =  5  l'r,i     -     °  y'7,10  =  2 
y'8,6  =  6   y,8,7=  7  y'8,10  =  3  y,9,10  =  6 

y'io,i =  20 

6.  Note  that  the  flow  in  all  sink  arcs  and  source  arcs  is  at 
arc  capacity,  therefore,  a  feasible  flow  exists.  Dropping 
nodes  0  and  11  plus  the  associated  arcs,  maximize  flow 
from  node  1  to  node  10  using  the  Ford  Fulkerson  algorithm. 
In  this  case  the  feasible  flow  computed  is  the  maximum 
flow.  Then,  using  the  change  of  variables  equations  the 
resulting  maximal  flow  on  G  is  determined.  Finally  G 
is  simplified  by  eliminating  arcs  ag  g  and  ag  q.  The 
maximal  flow  is  20.  The  flow  in  each  arc  is  as  follows: 

*2,4  =  8  ^2,6  =  4  *3,4  =  3  *3,5  =  5 

*4,5  =  7  y4,6  =  7  *5,4  =  3  *5,7  =  7 

*>,10  =  5  ^6,4  =  °  ^6,5  =  2  ^6,7  =  3 

*6,10  =  9  y7,5  =  ]  *7,6  =  3  ^7,10  =  6 

The  flow  in  node  5  is  9  and  in  node  6  is  9. 

Computer  Program  Description 

Since  it  is  often  the  case  that  it  is  necessary  to  solve  the 
maximal  flow  in  a  much  larger  and  more  complicated  network,  a  computer 
program  for  the  solution  of  the  generalized  maximal  flow  problem  as 
presented  in  this  thesis  has  been  written.  The  program  is  written  in 
FORTRAN  IV  for  use  on  an  IBM  360  computer. 
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START 

— EZ 

READ  IN 
NETWORK 
G 


D 


ALL  ARCS  INCIDENT 

TO  AND  INCIDENT 

FROM  THE  SAME  SET 

OF  NODES  ARE 

COMBINED 


NUMERICAL  VALUE 
OF  THE  TERMINAL 
NODE  IS  INCREASED 
TO  PERMIT  EXPAN- 
SION OF  THE  NETWORK 


COMBINE  MULTIPLE 
ARCS  JOINING  TWO 
NODES  IN  OPPOSITE 
DIRECTIONS 


TRANSLATE  BOUNDED 

NODES  INTO 

BOUNDED  ARCS 


COMPUTE  CAPACITIES 

FOR  NEW  SOURCE  AND 

SINK  NODES  IN  G' 


I 


CREATE  G' 


I 


GENERAL  FLOWCHART 
30 


*-<- 


¥. 


CREATE  CHAIN 

FROM  SOURCE  TO 

SINK 


LOAD  PATH 
ARRAY 


I 


SINK 
-L±±<         LABELED 
? 


COMPUTE  MAXIMUM 
FLOW  INCREASE 


FLOW  IS  INCREASED 
ALONG  PATH 


YES 


ELIMINATE  ALL 

ARCS  FROM  SOURCE 
AND  SINK  NODES 
OF  G' 


NO 


MO 


TRANSFORM  FLOW 
FROM  G'  TO  G 


m 


TRANSFORM  BOUNDED 

ARCS  TO  BOUNDED 

NODES 


I 


SEPARATE  APPRO 

PRIATE  ARCS  AND| 

FLOWS  TO  TWO 

OPPOSITE  ARCS 

AND  FLOWS 


I 


CHANGE  VALUE  OF 
TERMINAL  NODE  AND! 
COMPUTE  MAXIMUM 
FLOW  FOR  G' 


I 


RINT  ARC  FLOWS 
AND  MAXIMUM 
FLOW  FOR  G 


GENERAL  FLOWCHART  CONTINUED 
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The  maximal  flow  in  generalized  networks  1s  computed  using 
the  methodology  outlined  in  this  chapter.  In  addition,  a  set  of 
routines  is  added  to  the  program  which  will  permit  the  user  to 
transfer  several  networks  to  input  data  directly  without  first 
reducing  the  networks  to  basic  networks  by  hand  or  auxiliary  methods. 
However,  as  the  program  is  written  each  network  must  have  a  single 
source  and  sink  and  all  undirected  arcs  are  to  be  represented  by 
two  oppositely  directed  arcs.  The  program  has  been  checked  out  and 
functions  properly. 

The  present  program  is  written  for  a  50  node  network.  In 
order  to  solve  larger  networks  the  size  of  the  arrays  in  the  dimension 
statement  must  be  increased.  The  general  flow  chart  of  the  entire 
program  is  given  on  pages  30  and  31.  For  the  interested  reader  a 
detailed  flowchart  is  presented  in  Appendix  A.  A  Fortran  listing 
of  the  program  is  contained  in  Appendix  B  along  with  the  solution 
to  the  previous  example. 

Conclusions  and  Recommendations 

The  general  maximal  flow  problem  can  be  solved  more  efficiently 
by  use  of  the  method  developed  by  H.  Greenberg.  The  advantages  of  the 
technique  are  evident  when  solving  small  problems  by  hand.  The  tech- 
niques presented  provide  a  systematic  approach  to  obtaining  a  solution. 
The  problem  of  computing  the  maximal  flow  in  a  mixed  network  in  which 
the  flow  in  the  undirected  arcs  are  bounded  below  by  a  positive 
number  cannot  be  solved  by  this  method.  In  fact,  to  the  best  of  my 
knowledge  this  problem  is  unsolved  in  the  general  case  except  by 
complete  enumeration  techniques. 


32 


BIBLIOGRAPHY 


1.  Claude  Berge,  The  Theory  of  Graphs  and  Its  Applications,  John 
Wiley  and  Sons,  1958. 

2.  C.  Berge  and  A.  Ghouils-Houri ,  Programming,  Games  and  Transporta- 
tion Networks,  John  Wiley  and  Sons,  1965. 

3.  R.  G.  Busaker,  T.  H.  Saaty,  Finite  Graphs  and  Networks,  McGraw- 
Hill,  1965. 

4.  G.  B.  Dantzig,  Linear  Programming  and  Extensions,  Princeton 
University  Press,  1963. 

5.  G.  B.  Dantzig,  D.  R.  Fulkerson,  On  Max  Flow  Mi n  Cut  Theorem 
of  Networks,  Rand  Research  Memorandum,  R  M  -  1418  -  1 ,  1955. 

6.  G.  B.  Dantzig,  D.  R.  Fulkerson,  Computation  of  Maximal  Flow 
in  Networks,  Rand  Research  Memorandum,  R  M  -  1489,  1955. 

7.  L.  R.  Ford  and  D.  R.  Fulkerson,  Flows  in  Networks,  Princeton 
University  Press,  1962. 

8.  D.  R.  Fulkerson,  A  Network  -  Flow  Feasibility  Theorem  and 
Combinatorial  Applications,  Canadian  Journal  of  Mathematics, 
Vol.  II  (1959),  pp.  440-451. 

9.  Gale,  David,  A  Theorem  on  Flows  in  Networks,  Pacific  Journal  of 
Mathematics,  Vol.  7  (1957),  pp.  1073-1082. 

10.  Greenberg,  H.,  On  Generalized  Maximal  Flow  in  Networks,  unpublished, 

11.  Ford,  L.  R.  and  Fulkerson,  D.  R.,  Network  Flow  and  Systems  of 
Representatives,  Canadian  Journal  of  Mathematics,  Vol.  10  (1958), 
pp.  78-85. 

12.  _^ ,  A  Simple  Algorithm  for  Finding  Maximal  Network  Flows 

ancl  Applications  to  the  Hitchcock  Problem,  Canadian  Journal  of 
Mathematics,  Vol.  9  (1957),  pp.  210-218. 


33 


APPENDIX  A 

DETAILED  FLOW  CHARTS 

This  appendix  contains  the  detailed  flow  chart  of  the  computer 
program  in  addition  to  a  description  of  variables  used.  The  statements 
are  written  in  FORTRAN  IV  with  arrows  indicating  the  logical  flow 
through  the  system.  The  input  and  output  statements  are  not  enclosed. 
Description  of  Variables 

Variable  Name  Description 

C(I,J)      Capacity  of  the  directed  arc  from  node  I  to  node  J 

B( I , J)      Lower  bound  of  the  directed  arc  from  I  to  node  J. 

Y(I,J)      Initially  used  to  read  in  data. 

Y(I,1)      The  number  of  the  node  the  arc  a.  .  is  incident 
from.  *J 

Y(I,2)      The  number  of  the  node  the  arc  a.  .  is  incident  to. 

Y(I,3)      Upper  bound  of  arc  a.  .. 

Y( I ,4)      Lower  bound  of  arc  a.  .. 

Y( I , J)  The  value  of  the  flow  in  a  directed  arc  from  I 
to  J  after  input  data  is  transmitted  to  arrays 
C(I,J)  and  B(I,J). 

XL(I)       Dummy  array  used  to  compute  values  of  the 

capacities  of  the  source  and  sink  arcs  for  G' 
and  then  the  label  of  the  nodes  in  the  label 
routine. 

NL(I )       Dummy  array  used  in  the  label  routine  as  a  check 
to  determine  if  all  nodes  have  been  checked  and 
then  to  store  the  path  used  to  increase  flow. 

XNC(I.l)  Number  of  the  node  that  is  bounded. 

XNC(I,2)  Upper  bound  on  a  flow  through  a  node. 

XNC(I,3)  Lower  bound  on  flow  through  a  node. 

XNC(I,4)  Value  of  flow  through  a  node. 

NN  Number  of  networks. 
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N         Number  of  nodes  in  the  original  problem. 

M         Value  of  the  designated  number  of  the  sink  node. 

KX        Value  of  the  designated  number  of  the  source  node 

NA         Number  of  arcs  in  the  original  problem. 

NC         Number  of  bounded  nodes. 

Input  Card  Format 

All  entries  on  the  card  format  in  this  section  are  right, 
justified  and  are  fixed  point  entries  if  the  variable  name  begins 

with  I,  J,  K,  L,  M,  or  N  and  are  floating  point  entries  for  all  other 
variable  names. 
Card  #1 


Column  Name 

1-4  NN 

Card  #2 

1-4  N 

4-8  NA 

8-12  NC 

Card  #3  through  Card  NA  +  1 

1  -  10  Y(I,1) 

11  -  20  Y(I,2) 

21  -  30  Y(I,3) 

31  -  40  Y(I,4) 
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Card  #NA  +  3  through  Card  NC  +  3 
Column       Name 


1 

-   10 

XMC(I.l) 

11 

-  20 

XNC(I,2) 

21 

-  30 

XNC(I,3) 

Detailed  Flowchart 

The  detailed  flowchart  of  the  program  is  given  in  the  next 
11  pages. 
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KX  =  1 

M   =  N  +  NC  +  ; 

RET  =  0 


I 


DO  I     = 

Ml     /  T  \      _ 

! 
1,50 

ViL\  I  )     -    \J 

XL(I)   =  0                j 

ir 

DO  J       =  1,50 

B(I,J)  ■ 
C(X,J)  = 

U 
0 

DO  X       m   1,NA 


L  =  Y(I,1)  +  1 
K  =  Y(I,2)  +  1 


I 


C(L,K)  =  Y(I, 3} 
+  C(L,K) 


B(L,K)  =  Y(T,4) 

+  B(L,K) 


DETAILED  FLOWCHART 
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L  =  N  + 

1 

K  =  N  +  NC  +  1 

v 

DO  I  =  1,M 

v 

C(I,K)   =  C(I,L) 
B(I,K)  =  B(I,L) 

i 

' 

C(I.L)  =  0 

dU»l;  - 

u 

C(J,I)  =  -  C(J,I) 


DO  I  =  1,M 
DO  J  =  1  ,M 


NO 


C(I,J)  =  -C(I,J) 


DETAILED  FLOWCHART 
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/ 


< 


/     NC  >  0 


L  =  N  +  1 


DO  I  =  1,NC 
-M  K  =  XNC(I,1)  +  1 


DO  J  =  1  ,M 


NO 


NO 


s  \ 

^'C(K,J)  >  0     X. 


C(K,J)  =  -  C(K,J) 


C(K,L)  =  XNC(I,2) 
B(K,L)  =  XNC(I,3) 
L  =  L  +  1 


DETAILED  FLOWCHART 


YES 


C(L,J)   =  C(K,J) 
B(L,J)   =  B(K,J) 
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C(M-1,2)  =  .9999 
MX  =  M-l 


I 


DO  I  =  2,MX 


C(I,M)  =  -XL(I) 


NO 


DO  I  =  1,50 
XL(I)=  0 


YES 


C(l.I)  =  XL(I) 


DETAILED  FLOWCHART 
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NL(I)  =  0 


>   DO  J  =  KX,M 


-A_ 


XL(KX)  =  KX 


DO  L  =  KX,M 
I  =  L 


NO 


XL(J)  '  I 


NL(J)  =  XL(J) 


I 


LN  =  L 


DO  K  =  1  ,LN 


I  =  K 


NO                   NO 
VK 


XL(J)  =  I 


DETAILED  FLOWCHART 
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NO 


I  =  1 
NL(I)  =  M 

K  =  NL(I) 


1  =  1  +  1 

NL(I)  =  XL(K) 

K  =  NL(I) 


v   X  =  9999.9 


DO  I  =  1,M 
K  =  NL(I) 
L  =  NL(I+1) 


X  =  C(L,K)-Y(L,K) 


X  =  Y(K,L) 


DETAILED  FLOWCHART 
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A      <" 


(    r    ^ 

DO  I  =  1,M 
K  =  NL(I) 
L  =  NL(I) 

-/ 

vl 

^NO 

ZERO  OUT  ARRAYS 
NL  AND  XL 

C          L  >  °          > 

JYES 

NO 

Y(L,K)   = 
Y(L,K)   +  X 

^C(K,L)  >o\ 

YES 

Y(K,L)  = 
Y(K,L)   -  X 

1 

:: 

d  u m 


c 


STOP 


NO 


MJ  =  M  -   1 


KX  =  2 
M  =  M  -   1 
C(M,KX)   =  0 


DETAILED  FLOWCHART 
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^v_. 


NO 


DO  I  =  KX,M 
DO  J  =  KX,M 


C(I,J)  =  C(I,J) 
+  B(I,J) 

Y(I,J)  =  Y(I,J) 
+  B(I,J) 


j   C(K,L)  =  0 
j   Y(K,L)  =  0 


YES 


YES 


■*n 


DO  I  =  1,NC 

K  =  XNC(I,4)+1 


XNC(I,4)  =  Y(K,L) 


C(K,J) 
Y(K,J) 
Y(L,J) 


i  C(L,J)  = 


C(L,J) 
Y(L,J) 
0 
0 


DETAILED  FLOWCHART 
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DO  I  =  KX.M 
DO  J  =  KX,M 


C(J,I)  =  -  C(J,I) 

C(I,J)  =  C(I,J)  + 
B(J,I) 


I 


B(I,J)  = 
B(I,J)  +  C(J,I) 


Y(J,I)  =  B(I,J) 
-  Y(I,J) 


L  =  N  +  1 

K  =  N  +  NC  +  1 


DO  I  =  KX,M 
C(I,L)=  C(I,K) 
B(I,L)=  B(I,K) 


I 


Y(I,L)  =  Y(I,K) 
C(I,K)  =  0 


DETAILED  FLOWCHART 
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Y(J,I)  =  B(J,I) 

Y(I,J)  =  Y(I,J) 

+  B(J,I) 


Y(I,J)  =  B(I,J) 


PRINT 

SOLUTION 

OF  NETWORK 

FLOW 


STOP 


DETAILED  FLOWCHART 
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APPENDIX  B 

FORTRAN  -  IV  LISTING 
This  appendix  contains  the  complete  listing  of  the  computer 
program.  The  problem  illustrated  in  Chapter  IV  was  used  to  display 
the  output  format.  Comment  cards  are  included  for  clarification  and 
identification  of  the  main  sections  of  the  program. 
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DIMENSION  C.(  50,50)  ,B(  50,50  )  ,Y(  50,  50),  XL  (50)  ,NL(  50), 
1XNC( 50,4) 
READ(5,801)  NN 
801  FORMAT* 14) 

DO  800  IU=1 ,NN 
READ(5tl)  N,NA  , NC 

1  F0RMATI3I4) 
00  21=1, NA 
READ(5.3)(  Yd  ,J),  J=l,4) 

3  F0RMAT(4F10.4  ) 

2  CONTINUE 

DO  8  I=1,NC 

READ (5,4) (XNC( I, J) ,J=1,3) 

4  F0RMAT(3F10.4) 
8  CONTINUE 

KX  =  1 

M=  N  +  NC  +  2 

RET  =  0 

DO  5  1=1,50 

NL(I)=0 

XL(I)  =0 

DO  5  J=l,50 

B(  I, J)=0 

C(I,J)=0 

5  CONTINUE 

C  THE  CAPICITY  MATRIX  IS  LOADED  AND  ALL  ARCS  INCIDENT  TO 
C  AND  INCIDENT  FROM  THE  SAME  NODES  APE  COMBINED 

DO  6  I=1.NA 

L=  Y(I,1)  +  1. 

K=  Yd, 2)  +  1. 

C(L,K)=  Yd, 3)+  C(L,K) 

B(L,K)=  Yd  ,4)4-  B(L,K) 

6  CONTINUE 

C  THE  NUMBERICAL  VALUE  OF  THE  TERMINAL  NODE  IS  CHANGED  TO 
C  ALLOW  FOR  EXPANSION  OF  THE  NETWORK 

L=N+1 

K=N+NC+1 

DO  11  I=1,M 

C( I,K)=C(I,L) 

8( I,K)=B(I,L) 

C( I,L)=0. 
11  B( I,L)=0. 
C  COMBINE  MULTIPLE  ARCS  JOINING   TWO   NODES  IN  OPPOSITE 
C  DIRECTIONS. 

DO  15  1  =  1, M 

DO  15  J=1,M 

IF(C(I ,J) .GT.  0..AND.  C(J,I).GT.  0.)  GO  TO  16 

GO  TO  15 

16  IF(C(I,J) .LE.B( J, I  )  )  GO  TO  17 
C( I, J)=  C(I , J)-B( J, I) 

B(  I,J)=  Bd,J)-C(  J, I) 
C( J, I )=-C( J, I) 
GO  TO  15 

17  C(J,I)=C(J,I)-  B(I,J) 
B( J,I)=B( J, I)-  CII.JI 
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C( I, J)=-CU ,J) 

15  CONTINUE 
C  TRANSLATE   BOUNDED  NODES  INTO  BOUNDED  ARCS 
IF(NC)155,155,156 
156  L  =  N+1 

DO  34  1=1, NC 
K=  XNC( I  ,I)  +  1. 
DO  32  J=l,  M 
IF(C(K,J))  32,32,33 

33  C(L,J)=  C(K,J) 
B(L,J)=  B<K,J> 

IF(C(K,J)  .GT.O.)  C(K,J)=-C(K,J) 
32  CONTINUE 

C(K,L)=  XNC( 1,2) 
8(K,L)=XNC< 1,3) 
L  =  L+1 

34  CONTINUE 
155  CONTINUE 

C  THE  FOLLOWING  OPERATION  COMPUTES  THE  CAPACITIES  OF  THE 
C  ARCS  FROM  AND  TO  THE  SOURCE  AND  SINK  NODES  IN  G' 

DO  9  I  =  1,M 

SI=0 

SO  =  0 

DO  10  J=1,M 

Y( I, J)=0 

IF(C(I,J).GT.  0.)  C(I,J)=  C( I,J)-B( I, J) 

IF(C(If JJ.GT.O. )  SO=SO*B(I,J) 

IF(C(J,I).GT.O)  SI=SI+B(J,I) 
10  CONTINUE 

XL(  I)=SI-SO 
<5  CONTINUE 
C  GMS  COMPLETED  BY  ADDING  NEW  SINK  AND  SOURCE  ARCS 

C(M-1,2)=9999 

MX=  M-l 

DO  12  1=2,  MX 

IF(  XL(D)  13,12,14 

13  C( I,M)=  -XL(I) 
GO  TO  12 

14  C(1,I)=  XL( I) 
12  CONTINUE 

DO  31  1=1,50 
31  XL(I)=0 
20  XL(KX)=KX 
C  THE  LABEL  ROUTINUE  AS   PRESENTED   BY  FORD  AND  FULKERSON 
C  IS  USED  TO  DETERMINE  THE  FLOW  CAN  BE  INCREASED 

DO  210  L=KX,M 

I  =  L 

IF(XL( I) )93,93,96 
96  NL(I )=0 

GO  TO  22 
93  LN=L 

DO  94  K=1,LN 

IF(NL(K) )94,94,95 
95  I=K 

GO  TO  96 
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94  CONTINUF 
GO  TO  210 

22  00  21  J=  KX,  M 
IF(XLU) )23,23,21 

23  IF(C(I,J))  24,24.25 

25  IF(  Y(I,J).LT.  ClIiJII  XL(J)=I 
NL< J)=XL( J) 
GO  TO  21 

24  IF(C(J,I ))21,?1,26 

?6  IF(Y( J, I ) .GT.  0.)XL(J)=I 

NL( J)=XL( J) 
21  CONTINUE 
GO  TO  93 
210  CONTINUE 

IF(XL(M)  )27,  27,28 
C  A  PATH  FROM    SOURCE  TO   SINK   IS   DETERMINED  DIRECTLY 
C  FROM  THE  LABELLED  NODES 
28  1  =  1 

NL(I )=M 
K=NL(I) 
30  1=1+1 

NL(I )=XL(K) 
K=NL(I) 

IF(  K  .EQ.  KX)GO  TO  40 
GO  TO  30 
4C  X=9999.9 
C  THE  MAXIMUM   AMOUNT  THAT  THE  FLOW  CAN  RE  INCREASED 
C  ALONG  THE  PATH  IS  COMPUTED 
45  DO  41  1=1, M 
K=NL(I) 
L=NL(I+1) 
IF(L)50,50,  42 

42  IF(C(K,L))43,43,44 

43  IF((C(L,K)-  Y(L,K)).LT.  X)  X=  C( L, K )-Y ( L, K  ) 
GO  TO  41 

44  IF{  Y(K,L).LT.  X)  X=  Y(K,L) 
41  CONTINUE 

50  CONTINUE 

C  THE  FLCW  IS  INCREASED 
DO  51  1=1, M 
K=NLU  ) 
L=NL(I+1) 
IF(L)60,60,  52 

52  IF(C(K,L)I  53,53,54 

53  Y(L,K)=  Y(L,K)+  X 
GO  TO  51 

54  Y(K,L)=  Y(K,L)  -X 

51  CONTINUE 

60  CONTINUE 

DO  61  1=  KX,M 
NL(I )=0 

61  XL(I)=0 
GO  TO  20 

C  THE  SINK  CANNOT  BE  LABELLED  AND  A  CHECK  IS  MADE 
C  TO  DERMINE  IF  THE  FLOW  IS  FEASIBLE 
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27  IF(M.EQ. (N+NC+1) )  GO  TO  80 
MJ=M-1 

DO  70  1=  2,  MJ 

IF((C(1,I»-  Yd, I))  .GT.  .0001  .OR.  (C(I,M|- 
1Y(I,M)).GT.  .0001)  GO  TO  90 
7C  CONTINUE 
KX  =  2 
M=M-1 
C(M,KX)=0 
GO  TO  20 
C  THE  ORIGINAL  NETWORK  IS  RECONSTRUCTED   AND  FLOW  IS 
C  TRANSFERED  FROM  G«  TO  G 

80  CONTINUE 

DO  81  I=KX,M 

DO  81  J=KX,M 

IF(C(I ,J).LT.O.)  GO  TO  81 

C(  I,J)s  C(I ,J)  +  B( I.J) 

Y(  l,J)=Yd,J)+B<  It  J) 

81  CONTINUE 
IF(NC)157,157,158 

158  L=N+1 

DO  100  1=1. NC 

K=  XNCd,l)d. 

XNC< 1,4)  =  Y(K,L) 

C(K,L)=0. 

Y(K,L)=0 

DO  101  J=KX,M 

IF(C(L,J).LE.O)  GO  TO  101 

C(K,J)=C(L,J) 

Y( K, J)=Y(L,J) 

Y(L,J)=0 

101  C(L,J)=0 
L=L*1 

100  CONTINUE 
157  CONTINUE 

DO  102  1=  K.X.M 

DO  102  J=  KX.M 

IF(C(J,I )) 103,102,102 

103  C< J,I)=-C< J, I) 

C( I,J)=  Cd,J)*B( J, I) 

B(  I.J)*  Bd,J)+C< J, I) 

IF((B(I,J)-B(  J, I))  .GT.  Yd, J))  GO  TO  104 

Y( J,  I)=B(J,I ) 

Y(  I«  J)=Yd,J)+B(  J, I) 

GO  TO   102 

104  Y( J,I)=B(I ,J)-Y( I, J) 
Y(  I,J)=B(I,J) 

102  CONTINUE 
L=Nd 
K=N+NC«-1 

DO  105  I=KX,M 
C(  I,L)=C(I,K) 
B(  I,L)=B(I,K) 
Y(  I,L)=Y(I,K) 

105  C( I,K)=0. 
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WRiTE(6,86) 

86  FORMATflHl,  25X,  •  ARC  ■  ,  9X  ,  •  LOWFR  ROUND' , 12X , '^L OW • , 
114X, 'UPPER  SOUND' ,//) 

DO  82  1=2, M 
RET=RET+Y< I ,M) 
DO  82  J=2,M 
IF(C( I,J))82,82,85 
85  L=I-1 
K=J-1 
WRITE (6, 83)  L,K,B(I,J),Y(I,J),C(I,J) 

83  F0RMAT(24X,I2, IX, 12,  5X ,F1 0. 3 , 1 IX , F10 . 3, 10X , F 10. ^ , / ) 
82  CONTINUE 

WRITE(6,87) 

87  F0RMAT(1H0,24X, 'NODF' ,9X, 'LOWER  BOUND' , 1 2X, • FLOW • , 
114X, 'UPPER  BOUND' ,//) 

IF(NC) 161,161,160 

160  DO  89  1=1, NC 
K=XNC(I,1) 

WRITE<6,88)     K     , XNC ( I , 3 ) ,XNC ( I , 4 ) , XNC <  1 , 2 ) 

88  FORMAT(2  5X,I2,6X,F10.3,UX,F10.3,10X,F10.3,//) 

89  CONTINUE 

161  CONTINUE 
WRITE(6,84)   RET 

84  FORMAT(  25X,'MAX  FLOW  =',F10.4) 
GO  TO  91 

90  WRITE  (6,92) 

92  F0RMAT(25X, 'NETWORK  HAS  NO  FEASIBLE  FLOW') 

91  CONTINUE 

DO  117  1=1,50 

DO  117  J=l,50 

Y( I,J)=0 
117  CONTINUE 
80C  CONTINUE 

ENO 
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A  «  C 


LOWER  ROUND 


FLOW 


MPPFD  pni!MD 


1 

2 

1 

3 

L 

4 

2 

ft 

5 

A 
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5 

4 

5 

h- 

ft 

5 
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5 

6 

5 

7 

s 

A, 

o 

U 

6 

5 

6 

7 

6 

8 

7 

c; 

7 

ft 

7 

8 

NO! 

)E 

5 
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0.0 

8.000 

0.0 

3.000 

5.000 

3.000 

1.000 

3.000 

0.0 

3.000 

2.000 

0.0 

0.0 

1.000 

3.000 

1.000 

3.000 
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LOWFP  BOUND 

4.000 

4.000 


12.000 
8.000 
8.000 
4.000 
3.000 
5.000 
7.000 
7.000 
3.000 
0.0 
7.000 
5.000 
0.0 
2.000 
3.000 
9.000 
1.000 
3.000 
6.000 

FLOW 

9.000 

9.000 


Q99.000 

99O.000 

13.000 

4.000 

7.000 
10.000 
13.000 

7.000 
13.000 

R.000 

7.000 

^.000 

2.000 

8.000 

3.000 

9.000 

4.000 

8.000 
10.000 

UPPFR  POUND 

15.000 

10.000 
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20.0000 
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